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0\ Abstract 

o 

We discuss how the different estimates of elliptic flow are influenced by flow fluctuations and 
04 nonflow effects. It is explained why the event-plane method yields estimates between the two- 
^ ^ particle correlation methods and the multiparticle congelation methods. It is argued that nonflow 
QJ effects and fluctuations cannot be disentangled without other assumptions. However, we provide 
C/^ equations where, with reasonable assumptions about fluctuations and nonflow, all measured val- 
QQ ues of elliptic flow converge to a unique mean V2,pp elliptic flow in the participant plane. Thus, 
I the 20% spread in observed elliptic flow measurements from different analysis methods is no 
longer mysterious. 

D 

I ElUptic flow has proved to be very valuable for understanding relativistic nuclear colli- 

Q sions im m. However, different analysis methods give results which spread over a range of 

^ 20% 1(3]. A higher accuracy is now needed because when comparing to relativistic viscous hy- 

, ^ , drodynamic calculations, an uncertainty of 30% in the elliptic flow parameter V2 leads to an 

uncertainty of 100% in the ratio of shear viscosity to entropy [4 |. 
CSJ For simplicity we will write v{ ) instead of v„{ } and cos(...) instead of cos[n(...)], where n is 

the harmonic number of the anisotropic flow. The final equations are independent of n. 

Flow Methods: The two-particle cumulant method v{2) corTelates each particle with every 
other particle, and is defined as 



^ y{2) = V<cos(0i - 02)> , (1) 

where { ) indicates an average over all particles in all events. The four-particle cumulant method 
f-*) v{4) is defined as 

. ^ v{4) = (2<cos(0, - cf>2)f - <cos(0i + <^2 - 03 - 04)))'^^ . (2) 

X 

^ The Lee- Yang Zeros method v{LYZ) is also a multiparticle correlation. The event-plane estimate 

of anisotropic flow is defined as 

v{EP) = <cos(0 - ^'r))/R , (3) 

where the particle of interest is always subtracted to avoid autocorTelations. R is the event plane 
resolution correction which is determined from the correlation between the event plane vectors 



'Condensed from arXiv:0904.2315 [nucl-ex] 
Preprint submitted to Nuclear Physics A 



September 18, 2009 



of two independent "subevents" A and B. Methods of choosing the subevents are randomly, 
according to pseudorapidity or charge, or combinations of these. 

In the special case where the event plane comes from only one subevent the resolution cor- 
rection is the subevent resolution 

R = VWI^^^ ■ (4) 

The corresponding estimate of anisotropic flow will be denoted by v{subEP), or, more particu- 
larly, v{etaSub) or yjranSub), depending on how the events were divided. 

In the more general case when the event plane comes from the full event, one first estimates 
the resolution parameter Xs of the subevents by solving numerically the equation 

^Or.)= V(cos(Ta-Tb)), (5) 
where the function K is defined by |l5]|6l 

where /q and /i are modified Bessel functions. Generally, the resolution parameter is related to 
the flow through ;^'s = v '\/N/2. One then estimates the resolution parameter ;^f of the full event as 
X = Xs The resolution correction for the full event R is defined by R = "Rix) - f^iXs V2)- 
Fluctuations: Elliptic flow is driven by the initial eccentricity of the overlap almond 
This eccentricity fluctuates from one event to the other It is fluctuations which make (v) in 
the participant plane larger than in the reaction plane. The magnitude of flow fluctuations is 
characterized by cr,,, defined by 

(tI ^ {v^) - {vf, (7) 

where v is the flow in the participant plane vpp in the case of fluctuations in the participant plane. 
Flow methods involve various functions of v, which are also affected by fluctuations. 

We derive the effect of fluctuations on the various flow estimates, to order cr^. Using the 
definitions of v{2} and v{4} from Eqs. ([T]i and (|2]l, 

v{2)2 = (v2) = (v)2 + £r2 (8) 

and 

v{4f = (2(vf -(v^))''^<v>2-c.^ (9) 

Fluctuations increase v{2) and decrease v{4) compared to vpp. 

Nonflow Effects: Now we discuss nonflow effects while neglecting fluctuations. The two- 
particle azimuthal correlation gets contributions from flow and from other "nonflow" effects 

<cos(0i - 02)) = <v)' + 5 , (10) 

where 6 is the nonflow part. One expects that 6 varies with centrality like 1 /N, where is some 
measure of the multiplicity 15] O. 

Using Eqs. ([TJ and ^TU[ , one obtains, to leading order in 6 



v{2f^{vf+6. 
2 



(11) 



On the other hand, v{4) is insensitive to nonflow effects, and thus 



y{4) = <y> 



(12) 



Equations: We assume that to leading order in crj and 6, the contributions of nonflow and 



fluctuations are additive. Eqs. ([8]) and (111 yield 



Similarly, Eqs. ^ and ([12]) yield 



v{2f = {vf +6 + 0-1. 



(13) 



(14) 



Although this equation was derived for v{4} it should apply to all multiparticle values. As for the 
event-plane methods, one can show that 
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('0 + h)^ , 

where ;o,i is a shorthand notation for I() \(x\l2) and /o,i for /o,i 0(^^/2). The differences between 
the various V2{ ) estimates always scale like 6 + 2cr^. This shows explicitly that fluctuations and 
nonflow effects cannot be disentangled with only these measurements. Thus we have defined 
<^S + 2(rl. 

Application to Data: So far the equations have used generic fluctuation and nonflow param- 
eters. To apply the analytic equations to extract (y) in the participant plane from experimental 
data, we have assumed that the fluctuations in v have the same fractional width as the fluctuations 
of the participant eccentricity 

(Ty = {(Tj{e)) (v) . (17) 

A nucleon Monte-Carlo Glauber calculation was used to calculate the fractional standard devia- 
tion of Spait f9l. For the nonflow contribution we have taken the value from proton-proton colli- 
sions and scaled it down by the number of participants. The value of 6pp was obtained by integrat- 
ing the minimum bias p + p curves of Ref. ifTOl . Fig. 1, and it was found that 6pp - 0.0145 ifTTI . 
Thus for nonflow as a function of centrality we assume 



6 = 6pp 2/N. 



part 



(18) 



knowing that in a p + p collision there are two participants. 

The published STAR data [3 12| for the various methods are shown in Fig. [T] The upper 
lines are from "two-particle" correlation methods, and the lower line is from a multiparticle cor- 
relation method. The lower line values for V2{LYZ) are thought to be in the reaction plane, if 
the fluctuations are Gaussian |fT3]| . The line for V2{etaSub) is somewhat low for peripheral col- 
lisions because the gap in pseudorapidity reduces short-range nonflow correlations. Particularly 
puzzling is why the V2{EP) line is lower than the other two-particle methods. 

Correcting to (v) in the participant plane was done by using Eq. ( 13 1 for V2{2), Eq. ( 14 1 for 
V2{LYZ), Eq. ( 15 1 for V2{EP), and Eq. ( 16 1 for V2{ranSub) and V2{etaSub). The results are shown 
in Fig.|2] Since v'2{etaSub) is less affected by nonflow, the value of 5pp used for it was multiplied 
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Figure 1 : (Color online) The values of V2 from var- 
ious analysis methods vs centrality. Both the upper 
lines (3J and the lower line [12] are STAR data. 



Figure 2: (Color online) The data from Fig.[T|cor- 
rected to (vo) in the participant plane. 



by 0.5. In Fig.|2]the convergence of the two-particle, full event plane, and multiparticle results 
to one locus in the participant plane is remarkable. Even the shape of the V2{etaSub) curve has 
changed to match the others with only one additional parameter. Previously we took the spread 
in the values in Fig.[T]as an estimate of the systematic uncertainty. 

Summary: We have shown how the various experimental measures of elliptic flow are af- 
fected by fluctuations and nonflow, and we derived analytic equations which are leading order in 
erf. and 6. We have transformed published data to the participant plane using reasonable assump- 
tions for fluctuations and nonflow. The convergence of the various experimental measurements 
is remarkable. The convergence of the methods essentially fixes the value of crtot from experi- 
mental data, but the separation into fluctuation and nonflow parts is not unique. To avoid both, 
better results for multiparticle correlations are needed. 
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